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The electron addition spectrum A'^{k,Lo) is obtained analytically for the one-dimensional (ID) 
supersymmetric t-J model with interaction. The result is obtained first for a small-sized system 
and its validity is checked against the numerical calculation. Then the general expression is found 
which is valid for arbitrary size of the system. The thermodynamic limit of (fc, oj) has a simple 
analytic form with contributions from one spinon, one holon and one antiholon all of which obey 
fractional statistics. The upper edge of A'^{k,ui) in the {k,Lo) plane includes a delta-function peak 
which reduces to that of the single-electron band in the low-density limit. 
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The concept of spinons and holons, both of which obey 
the fractional statistics has turned out to be use- 
ful in approaching to ID electron systems. In terms of 
these quasi-particles one can inquire into not only the 
low-energy and low- wavelength limit, but the global fea- 
ture of the dynamics. Hence special interest has been 
cherished in the global dynamics from both theoretical 
and experimental points of view. For example, angle 
resolved photoemission Q| has revealed some evidence 
of the spin-charge separation by resolution of holon and 
spinon contributions. On the theoretical side, numerical 
studies have been performed for the ID t-J model for a 
small number of lattice sites and some structures have 
been ascribed to spinons and holons. For deeper under- 
standing of the overall dynamics, demand is growing for 
analytic theory which can go to the thermodynamic limit. 
Partly analytic theory is available for the single-particle 
spectral functions of the t-J model in the J — > limit 

. A notable feature is that a satellite band is observed 
whose intensity is comparable to that of the main band. 
It is natural to ask how the finite J influences the dy- 
namics. 

In the supersymmetric t-J model with interac- 
tion Q , spinons and holons appear in the simplest man- 
ner. In fact exact thermodynamics for the model ^ can 
be interpreted in terms of free spinons and holons. Ha 
and Haldane |^ analyzed numerical results for dynam- 
ics in finite-sized systems, and found that only a few 
number of elementary excitations contribute to spectral 
functions. They proposed a momentum-frequency region 
where each spectral function takes nonzero values in the 
thermodynamic limit, but they did not obtain the spec- 
tral functions themselves. Recently, exact results have 
been derived for a particular component and for a 
particular momentum range of the spectral weight . 

In this paper we report on the analytical result of the 
electron addition spectrum for the t-J model at zero 
temperature. The electron addition spectral function 



is relevant to the angle resolved inverse photoemission 
spectroscopy. Our result constitutes the first analytical 
knowledge for dynamical quantities of lattice electrons 
with no restriction on the system size, the density and 
the momentum-frequency range. Although we cannot 
provide the formal proof for the exactness, the analytic 
results show complete agreement with numerical results 
for finite systems with various sizes. Hence our result in 
the thermodynamic limit is also expected to be exact. 
We consider the supersymmetric t-J model given by 



-t,j (cl^Cja + h.C. 



V, 



(1) 



where Cia- is the annihilation operator of an electron at 
site i with spin cr, rii is the number operator and Si is 
the spin operator. The projection operator V excludes 
double occupation at each site. The transfer and ex- 
change energies are given by ty = Jijf^ = where 
Dij = (N/n) sin (tt {i — j) /N) with even TV being the 
number of lattice sites, and the lattice constant is unity. 
The electron addition spectral function with the ground 
state |0) is defined by 

A+{k,u) =^|(i.;iVe + l|cL|0;A^c)|' 

x5 {uj ~ E,{N, + I) + Eo{N,) + C) , (2) 

where iVc is the total electron number, C, the chemical 
potential, cj,^ — A^"^^^ X); cL^''^', and \v) denotes an 
eigenstate of the Hamiltonian with energy E,j. 

We first give our main result and then its derivation. 
The addition spectrum includes the dispersion relations 
es{q) of spinons, th_{q) of holons, and Cal?) of antiholons. 
They are given in units of t by £3(9) = li^s — q), 
eh(g) = q{vc + q) and ea(g) = q{2vc - q)/2, where Ws = tt 
and Vc = 7r(l — n) with n = N^/N = 2k-p/'i^- 
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Analytical expression of A'^{k,uj) with < k < 2tt 
consists of the following components: 

A+ik,uj) ^ An{k,u) + Ai^{k,Lo) + Au{k,uj), (3) 

where 



AR{k,uj) 



An 



dqh 



dqs 



27r-4feF 



dqa 



xd{k - fcp - - (?h - 9a) 
x6{uj - esiqs) - eh(gh) - eal^a)) 

^, £s(gs)g°'^£h(gh)g"'^£a(ga)g°'^ 
(gh + ga/2)2 



(4) 



and Ai,{k,uj) = AR(27r — k,uj). In Eq.(||) we have 
gs = 1/2, (7h — 1/2 and = 2, which correspond to 
statistical parameters of spinons, holons, and antiholons, 
respectively [||j7|. Thus the matrix element in Eq.(U) 
can be interpreted in terms of the fractional statistics as 
in the case of correlation functions in related continuum 
systems [fo|. The third component is given by 



Au{k,u) 



'ea(fc-2fcF) 



eo(fc) 



S{UJ - Wau(fc)), 



(5) 



which contributes only in the region 2fcF < k < 2t: ^2kp. 
Here eo{k) = k{Ti — fc/2) describes the spectrum of non- 
interacting electrons, and ujs.\]{k) = es{kp) + ^^{k — 2k-p). 
The A-\j{k,uj) can be regarded as contribution from an- 
tiholons with fixed-energies of spinons and holons. In 
the dilute limit (fcp 0), tJau(fc) as well as ea(fc — 2fcF) 
tend to eo(fc)- Hence the coefficient in Eq.(||) becomes 
unity, and Au(fc, lo) tends to the spectral intensity of non- 
interacting electrons. The other contributions AR^L(fc, ^) 
can be neglected in this limit. 

Figure 1 shows the spectral edges by solid lines to- 
gether with the spectral intensities for a finite system as 
explained later. The threshold behavior in AY{{k,Lo) is 
derived as follows: As the frequency approaches an up- 
per edge given by iVsik) = es{k — kp) with kp < k < 2fcF, 
A-[{,{k,u)) diverges as [ajs(fc) —uj]^^^^. Along the connect- 
ing upper edge for 2fcF < k < 2tt — 2k-p, ^^(fc,^) has 
a stepwise discontinuity such as 9{uJiM{k) — uj). Here 
9{x) is a Heaviside step function. On the other hand, 
near a lower edge given by LO]^{k) = e\^{k -\- Zk-p — 27r) 
with 2tt — Zkp < k < 2tt — 2kp, AB,{k,uj) behaves 
as [uj — (jjh(fc)]'^^^- The feature is in contrast to the 
result of Ref.4 where the spectrum is enhanced near 
uj\i{k). Finally near another lower edge given by WaL(fc) = 
ea(fc — kp) with kp < k < 2Tr ~ 3kp, there is a stepwise 
discontinuity such as 9{uj — ujsj^{k)). We note that the 
asymptotic behavior of A'^ik^u) is fully consistent with 
the conformal field theory Q. 

Let us describe derivation of A+(fc, ui). We represent a 
state vector \^p) in the i- J model by 



\i.)=Y.^{x\x^)llSrl[h]\F). 



(6) 



where \F) is the fully up-polarized state, = q^Cij, 
and /i| = Cij. The set of coordinates {x^^x^) = x rep- 



resents x 



with n = A^h + Ns- Here x'' is for the i-th hole and 
Xj is for the j-th down-spin electron. In order to derive 
eigenfunctions of the system, it is convenient to use the 
Sutherland mod el |ll| with the SU(1,1) supersymmetry 
as an auxiliary [[12[|13[ . The merit of using the Suther- 
land model is that much more is known about properties 
of the eigenfunctions than those for the t-J model. The 
eigenfunctions with the coupling parameter /3 = 1 cor- 
respond to such part of the set {ij;{x^, x^)} that belongs 
to the Yangian highest weight states (YHWS) with the 
SU(2,1) supersymmetry iQ. In making the correspon- 
dence one restricts the coordinates x to integer variables. 
It is convenient to introduce the complex coordinate by 
Zj — exp{2TTixj / L) where L {— N) is the length of the 
ring-shaped system. The YHWS can be expressed by 
polynomials of Zj where the degree of each variables lies 
in the range [-N/2,N/2]. The YHWS form a subset of 
the Fock space and the other states can be generated by 
successive actions of Yangian generators on the YHWS. 

The ground-state wave function for the i- J model with 
iVh + 1 holes and Ns = Ne/2 down-spin electrons is given 

by [| 



*Gs = n 



z^ 

ie[0,n] iei, 



n 



n (^«-^.)' (7) 



i<je[0,n] 



where we have introduced the interval Is = [iVh -I- 1, A^h + 
Ns] for the suffices i and j. We have assumed that both 
Ns and iVh are odd so that the ground state is non- 
degenerate. In the following we always take /3 = 1. We 
fix a hole position at Xq = in v^qs which corresponds 
to adding an up-spin electron at this site. In terms of the 
ground-state wave function ^'gs for the system with TVh 
holes and Ng down-spin electrons, we can represent the 
resultant state by 



GS 



m=0 



(-l)™e™(z)*GS, 



(8) 



where em{z) = J2ii<i2< -<i^ei ' ' ' ^^e elemen- 

tary symmetric function of order to. Note that the 
interval [0,n] in Eq.(|^) is replaced by another interval 
I=[l,n] in Eq.(|). 

The spectrum of the Sutherland model Hsu is conve- 
niently analyzed with use of a similarity transformation 
generated by O = U.<jei (^^ " ^jf U.ei z-^^'^+^-^^Z^ 
We obtain H = O'^HsuO as follows: 
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2 n 



1=1 



(9) 



where di is called the Cherednik-Dunkl operator |15|. It 
is known that the set {di} can be diagonalized simulta- 
neously with real eigenvalue Xi for each di. The resul- 
tant eigenfunctions are polynomials E\{z; (3) of and 
are called the non-symmetric Jack polynomials . The 
polynomial E\[z;l3) has a triangular structure with re- 
spect to a certain ordering on the set of {A} [ p^|jT^ ]; if 
one expands E\{z](3) in terms of monomials T\iZ^ \ the 
expansion coefficient is zero unless A ^ i', where ^ de- 
scribes the ordering relation in the set. 

Since we are dealing with identical particles, eigen- 
functions should satisfy the following conditions of the 
SU(1,1) supersymmetry: 

(i) symmetric with respect to exchange of z^'s; 

(ii) anti-symmetric with respect to exchange of zf's. 

By taking linear combination of E\(z;(3), we can con- 
struct another polynomial K\{z;0) with the SU(1,1) su- 
persymmetry [p^p^. We specify the set of momenta 
as (Ah, A^) = (Af, . • • , Ah,,^, Al, • • • , A^^J with A^ > • • • > 
A^j^ , A| > ■ • • > A^v^ . In this way we can parameter- 
ize K\{z;(3) by using A = (Ah,A''). At the ground state 
with 2Ns + 1 electrons we have A = Aqs — (^gS'^^Gs) 
with A^s - {{Ns - l)/2, [N, - l)/2, • • • , (A^. - l)/2), and 
A^GS = (^s-l,^s-2,---,0). 

We define the inner product {f,g)o for complex func- 
tions f{z) and g{z) as the constant term in the Laurent 
expansion of f{z)*g{z). As is clear from the definition of 
the transformation O, K\{z; P)0{z) is orthogonal with 
respect to the above inner product. In order to derive the 
norm of Kx{z;(3), we generalize the procedure taken in 
Rcfs. The result of lengthy calculation is given 

by 



(10) 



where {ExO, E\0)o denotes the norm of the non- 
symmetric Jack polynomials. We refer to the literature 
for the explicit form of the latter norm since it re- 
quires many lengthy combinatorial quantities. The quan- 
tity px = p\pI is given by 

n ^T^' n 'vV^ ' (11) 

where we use the intervals /s and /h = [l,Ah] for the 
suffices i and j. 

In order to calculate we need to know the 

expansion coefficient cx which appears in 



,(2)^-03 = ^CAifA(^;/3)0(2). 



(12) 



A+(fc,cj) 



: Ah ^'(5(60 



(13) 



where ^Ex — Ex ~ Eq ~ the chemical potential ( 
for the finite A^ is given by C/'"''^ = — (3fi^ — 6n + 
4)/12 - (n - 2)/(2A^) -I- l/(3Af2), and the prime means 
that the summation over A is restricted by the momen- 
tum conservation |A| = |Ags| + "m- with |A| = J2j ^^'^ 
k ^ 2TT{m+{Ns+l)/2)/N. The norm (*gs, *Gs)o can be 
obtained by a procedure similar to that described above. 

From the triangular structure inherited to Kx, the 
contributions to the excited states are limited to the cases 
where — A^g + (1'*'% 0^="'*'=). Here means the se- 
quence ■ ■ ■ ,p with the number of p's being a. Therefore 
the spin part can be parameterized by the single vari- 
able Ag. We first calculate Kx by a brute force for small 
systems (up to A^h = 5 and Ag — 15) in the small mo- 
mentum region m < (A"s — l)/2. In this region, we have 
•^^ = ^GS + ^- The coefficient Cx is zero if the partition 
ly contains s — (2,2) where s = {i,j) denotes a square in 
the Young diagram |l5| . We find that the results for Kx 
are expressed in the following form: 



-a'{s) + l + l3'{l'{s)-l) 
a{s) + l + I3'l{s) 



(14) 



where /?' = /3/(/3 + 1) = 1/2, a{s) = Vt - j (arm- 
length), a'{s) = j — 1 (arm-colength), l{s) = u'^ — i (leg- 
length) with v'j the length of the column, and l'{s) — 
i — 1 (leg-colength). Thus, we can parameterize v as 
V — (Ah, l-'^a"!^ 0^''"'*'°) in this small momentum region. 
Namely A^{k,uj) is determined by the three parameters 
As, Ah and Aa, which are related directly with momenta 
of the excitations. We have checked numerically that 
the expression given by Eq.(^ can in fact be extended 
beyond the small momentum region. Note that cx is in- 
dependent of As which, however, enters A^{k, uj) through 
the norm {KxO,KxO)q. 

We thus obtain the finite-size version of Eq. ^ : 



A 



Ik{X)5{u: - AEx) 



(15) 



Using the coefficient cx, A+(fc,Ci;) can be expressed as 



and similar results with the suffix R replaced by L and 
U. In Eq.(p^ we have 

r ^ 1 r(Ah - l/2)r(Ah + Ah/2) 

^ 2 (r(i/2))2 r(Ah)r(Ah + (A^h + 1)/2) 

Aa(iVh - Aa + 1) 

''(2Ah + Aa-l)(2Ah + Aa-2) 
r(As + 1/2) r(A^/2 - As) 

r(As + i) r((A^ + i)/2-As)' 

with Nk/{2-K) = Ah + Aa + As + (A"s - l)/2, and 



(16) 



3 



27r 



2 r 



As(iV-l-2As 



+Ah(2Ah + - 2) + (Aa - 1){N^ - Aa) 



(17) 



The triangularity leads to < Ag + Ah < (iVg — 1)/2. The 
parameter Aa varies from zero to Ny^. For the case where 
As > (A^s ~ l)/2, we use the reflection symmetry about 
k = TT. We obtain /l(A) = /r(A) where the components 
in A = (As, Ah,Aa) are given by As = iVs - As, and by 
the relation A*^ + {0^'--^% Ah) = Xqs + (1^''). We 

obtain the allowed range < As + Ah < (iVg + l)/2. 

Finally the case where As = {Ng — l)/2 should be con- 
sidered separately. In this case, the triangular struc- 
ture requires that c\ should be unity, and that A*^ — 
^GS + (1^%0^''-^'') with 1 < Aa < A^h- We obtain for 
this special case 



^u(A) 



r[(Aa + 2)/2]r[(7Vh - Aa + 2)/2] 



r[(Aa + i)/2]r[(iVh-Aa + i)/2] 

r[(A^s + Aa)/2]r[(jVs + iVh - Aa + 2)/2] 
■ r[(7Vs + Aa + l)/2]r[(7Vs + ^h - Aa + 2)/2] 



(18) 



where 7Vfc/(27r) = Aa + N^. The excitation energy AE\ 
is given by Eq.([l7|) in the case (As,Ah,Aa) = {{Ng — 
l)/2,l,Aa). Thus in terms of /r, Jl and I\j we obtain 
the finite-size version of A'^{k,oj). 

In Fig. 1^, we present the result for N = 60, A'h — 29 
and Ns = 15. We have checked the validity of Eqs.(p^) 
and ( p^ ) by comparison with numerical results up to 
TV = 16 1 21 1 . In the special case = (0^''), momenta 
of the holon and the antiholon are both zero, and we ob- 
tain a form different from Eq.(|l6|). However, this case 
can be neglected in the thermodynamic limit. Following 
the same procedure as that in the spinlcss Sutherland 
model 1^,^, we obtain the expressions (^) and (||) in 
the thermodynamic limit. 

We hope that the quasi-particlc structure discussed in 
this paper is found in future experiment of inverse pho- 
toemission. 
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FIG. 1. The electron addition spectrum A~^{k,u)) in the 
case of N = 60, iVh = 29 and = 15 with Fermi momentum 
fcp — 7r/4. The intensity is proportional to the area of each 
oval. The solid lines are determined by dispersion relations of 
the elementary excitations in the thermodynamic limit. 
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